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Abstract 

We present the construction of a planar drift chamber with wires stretched between 
two arms of a c-shaped aluminium frame. The special shape of the frame allows to 
extend the momentum acceptance of the COSY-11 detection system towards lower 
momenta without suppressing the high momentum particles. The proposed design 
allows for construction of tracking detectors covering small angles with respect to 
the beam, which can be installed and removed without dismounting the beam-pipe. 
For a three-dimensional track reconstruction a computer code was developed using 
a simple algorithm of hit preselection. 
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1 Introduction 



The drift ctiamber wliicli is described in this report was built for the COSY-11 
experimental facility operating at the Cooler Synchrotron COSY-Jiihch [1]. 
The COSY-11, described in details in Ref. [2] and shown schematically in 
Fig. 1, is a magnetic spectrometer for measurements at small angles and is 
dedicated to studies of near-threshold meson production in proton-proton and 
proton-deuteron collisions (see e.g. Refs. [3,4,5]). It uses one of the regular 
COSY dipole magnets for momentum analysis of charged reaction products 
originating from interaction of the internal COSY beam particles with the 
nuclei of a cluster beam target [6]. Trajectories of positively-charged particles 
which are deflected in the dipole magnet towards the center of the COSY-ring 
are registered with a set of two planar drift chambers indicated as Dl and 
D2 in Fig. 1. These chambers cover only the upper range of momenta of the 
outgoing particles what suffices to measure e.g. two outgoing protons in the 
pp — > ppT]' process [5] or p — p — tracks in the pp ppK^K~ reaction 
close to threshold [4]. For tracking positively charged pions appearing in near- 
threshold reactions such as pp —>■ ppix^Ti^ with momenta by a factor of m.^lm.p 
smaller then the proton momenta it was necessary to extend the COSY-11 
momentum acceptance towards smaller values. 




Fig. 1. COSY-11 detection system. Indicated particle tracks correspond to a 
■pp — > ppvT+vr" event simulated for the beam momentum exceeding by 1 MeV/c 
the threshold momentum equal to 1219 MeV/c. The two outgoing protons are reg- 
istered with a set of two chambers — Dl and D2 and the 7r+ trajectory is measured 
with the presently introduced chamber — D3. The scintillation hodoscopes SI and 
S2 are used as start detectors and S3 hodoscope placed in a distance of 9.4 m 
from SI and not indicated in the figure is used as stop detector for time of flight 
measurements. 

Another important purpose was the detection of positively charged kaons prior 
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to their decay what is especially important for the measurement of the pp 
ppK^K~ close to threshold due to its small cross section on the level of a few 
nanobarns. For this, an additional drift chamber was built and installed in the 
free space along the COSY-11 dipole magnet (see chamber D3 in Fig. 1) . The 
main requirement for the chamber was a shape of the supporting frame which 
would not interfere with high momentum particles that are registered in the 
detectors Dl and D2. This demand was fulfilled by choosing the frame of a 
rectangular form but with one vertical side missing, called c-shaped frame since 
it resembles the character c. The main design characteristics of the chamber 
are given in section 2. It was also essential, that the chamber allows for a 
three-dimensional track reconstruction in order to determine the momentum 
vectors of registered particles at the target by tracing back their trajectories 
in the magnetic field of the COSY-11 dipole magnet to the nominal target 
position. Therefore three different wire orientations were chosen and a track 
reconstruction program was developed. This program is described in section 3. 
The chamber calibration and results of the track reconstruction are presented 
in section 4. 



2 The chamber construction 



The sensitive chamber volume is built up by hexagonal drift cells (see Fig. 2) 
identical with the structure used in the central drift chamber of the SAPHIR 
detector [7]. In this type of cells the drift field has approximately cylindrical 
symmetry, and thus the distance-drift time relation depends only weekly on 
the particles' angle of incidence. In order to minimize the multiple scattering 
on wires, gold-plated aluminium [8] is used for the 110 /xm-thick field wires, 
whereas the sense wires are made of 20 /im-thick gold-plated tungsten. 



field 
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sense 
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Fig. 2. Hexagonal cell. 



The cells are arranged in seven detection planes as indicated in Fig. 3 showing 
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one of two parallel aluminium endplates between which the wires are stretched. 
Three detection planes (1, 4 and 7) contain vertical wires, two planes (2, 3) 
have wires inclined at —10° and the remaining two planes (5, 6) contain wires 
inclined at +10°. This arrangement makes it possible to reconstruct particle 
trajectories in three dimensions, also in cases of multi-track events. The rel- 
atively small inclination of the wires was chosen since the measurement of 
particle trajectories in the horizontal plane needs to be much more accurate 
compared to the vertical plane. This is due to the fact that the particle mo- 
menta are determined on the basis of the deflection of their trajectories in an 
almost uniform vertical magnetic field in the 6 cm gap of the COSY-11 dipole 
magnet. 

Each endplate contains bore holes of 04 mm with inserted feedthroughs for 
mounting the wires. The drilling was done with a CNC jig boring machine 
assuring a positioning accuracy better than ±0.01 mm. The holes for the field 
wires form equilateral hexagons with a width of 18 mm. In the middle of these 
hexagons there are the openings for the sense wires. The resulting spacing of 
the sense wires in the planes with vertical wires is again equal to 18 mm and in 
the detection planes with wires inclined by ±10° is equal to 18 mmx cos 10° 
= 17.73 mm. The chamber contains 3 x 80 = 240 cells with vertical wires, 
2 X 76 = 152 cells with wires inclined by —10° and 2 x 77 = 154 cells with 
wires inclined by ±10°. 
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Fig. 3. Upper end-plate with indicated positions of openings for feedthrough used 
for mounting the wires. 



The disturbance of the electric field in the drift cells of a given plane by the 
cells of the neighboring planes was investigated using the Garfield code [9]. 
With the chosen spacing of the detection planes equal to 28.8 mm (see Fig. 
3) the resulting corrections to the distance-drift time relation determined as 
a function of position along the sense wires are smaller than 0.05 mm [10] 
for distances from the sense wire smaller than half of the cell width (9 mm). 
These corrections were neglected in view of the expected precision of track 
reconstruction in the order of 0.1 mm. Thus one can assume that the distance- 
drift time relation is the same for all the cells in a given detection plane, what 
simplifies the chamber calibration. 
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The two aluminium endplates for mounting the wires are 15 mm thick and are 
supported by two c-shaped frames made out of 20 mm thick aluminium (see 
a three-dimensional view in Fig. 4). The frames hold the total load of about 
2.4 kN originating from the mechanical tension of ~0.2 kN (= 546 x 0.3 N) 
of the sense wires and ~2.2 kN (= 2198 x 1 N) of the field wires. Prior to 
mounting the wires, the frames were pre-stressed with a force corresponding 
to the tension of wires using led bricks uniformly distributed on the upper 
endplate. This caused a deflection of the endplates at the free ends of 0.7 mm 
each. The ends of the pre-stressed endplates were fixed together using a steel 
plate and the bricks were taken away. After mounting the wires the steel 
plate was removed and we checked that the distance between the ends of the 
endplates did not change. 

Window frame U-shaped frames 




1500 mm 
1700 mm 



Fig. 4. Schematic three-dimensional view of the chamber frame for mounting the 
wires. The window for particles is 1500 mm wide and 390 mm high. 

For mounting the field wires brass feedthroughs with a 0150 fim inner opening 
are used. Since all field wires should have the same ground potential electrical 
contact of the feedthroughs with the aluminium endplates is assured by usage 
of a conductive glue. For the sense wires we used isolating feedthroughs made 
of Delrin, a polyacetal with a high dielectric strength of about 40 kV/mm, with 
inserted brass feedthroughs containing 050 fim openings for the positioning 
of the wires. A schematic drawing of the feedthroughs is shown in Fig. 5. The 
wires were strung between pairs of feedthroughs and were fixed by means of 
copper crimp pins. For sealing the feedthroughs and the crimp pins epoxy resin 
was used. The high voltage was connected to the sense wires through 1 MQ 
resistors soldered directly to the crimp pins on the lower endplate. The signals 
were read out through 1 nF coupling capacitors connected to the sense wires 
on the upper endplate (see Fig. 6). The capacitors and resistors were enclosed 
in hermetical volumes which were dried out with silica aerogel allowing to 
reduce leakage currents. 

For the window for particles penetrating the chamber 20 /^m thick Kapton foil 
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Fig. 5. Section of feedthrough for (a) field wires, and (b) sense wires. 



is used glued to an aluminium frame which is screwed on the support frame 
and sealed with 03 mm o-ring gasket. 
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Fig. 6. Section of the upper endplate. 



Each of the detection planes is equipped with five 16-channel preamplifier- 
discriminator cards [11] based on the Fujitsu MB13468 amplifier chip and the 
LeCroy MVL407S comparator chip. The cards are mounted directly on the 
chamber. The output pulses in the ECL standard are led by means of 30 m 
long twisted-pair cables to FASTBUS-TDCs of LeCroy 1879 type working in 
the common stop mode with the pulses from the discriminators as start signals 
and the trigger pulse as the common stop. 

As chamber gas, PIO mixture (90% argon and 10% methane) at atmospheric 
pressure is used. The gas flow through the chamber is about 12 1/h. The sense 
wire potential is -1-1800 V, whereas the field wires are all grounded. The gas 
amplification in the chambers is about 10^ and the discrimination threshold 
set in the preamplifier-discriminator cards corresponds to 3 • 10^ electrons. 
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3 Track reconstruction algorithm 

For the reconstruction of particle tracks a simple algorithm was developed and 
implemented as a computer code written in the C-language. The reconstruc- 
tion proceeds in three stages: 

• finding track candidates in two dimensions, independent for each orientation 
of the wires, 

• matching the two-dimensional solutions in three dimensions, 

• three-dimensional fitting in order to obtain optimal track parameters. 

In the first stage, track candidates are found in two-dimensions using hits in 
the detection planes with the same inclination of the wires. For this, all possible 
combinations of pairs of hits from two different detection planes are taken into 
account. Due to the left-right ambiguity of the track position with respect 
to the sense wire, for each pair of hits there are four straight-line solutions 
possible. These solutions are determined using an iterative procedure. First, 
the track distance to the sense wire is calculated from the drift time to drift 
distance relation d{t, 9) assuming the track entrance angle 9 as equal to 0°. 
Note that the time-to-distance function may depend on the angle 9 (see Fig. 
7) and this is taken into account when making a calibration. Thus, the track 
is defined by a pair of points lying in the corresponding sense wire planes and 
indicated as PI and P2 in Fig. 7. In the next step, the inclination of this track 
is taken for determining new values of the distances from the sense wires and 
constructing a new pair of points defining the track. These points are indicated 
as Tl and T2 in Fig. 7. This procedure can be repeated, however, we terminate 
it already after the second step since further steps give negligible corrections. 
For each track determined by this method, one subsequently finds hits which 
are consistent with it within a certain corridor along the track. In this way, hits 
from neighboring cells with respect to the selected two cells and from other 
detection planes with the same orientation of wires are taken into account in 
the reconstruction. For the width of the corridor we take ±1 mm which is a 
few times larger than the position resolution of a single detection plane. With 
this choice more than 90% of all hits are included in the reconstruction as it 
is discussed in the next chapter. 

In three dimensions a two-dimensional solution can be represented by a plane 
containing the track and parallel to the corresponding sense wires. Planes rep- 
resenting the two-dimensional solutions for three different directions of sense 
wires used in the chamber should intersect along a common line correspond- 
ing to the particle track. However, due to the limited spatial resolution of the 
chamber, the planes intersect along three different lines (see Fig. 8). Two lines 
— /i and I2 correspond to the intersections of the inclined planes with the 
vertical plane and the third line — Z3 is an intersection of the inclined planes. 
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Fig. 7. Illustration of iterative procedure of fitting a straight line to two hits. 

In order to quantify the consistency between the two-dimensional solutions, 
we determine the distances hi and /i2 between the crossing points of the lines 
h, I2 with the first and the seventh detection plane, respectively (see Fig. 8). 
For the three-dimensional reconstruction, we accept only the combinations of 
two-dimensional solutions for which the distances hi and /i2 arc both smaller 
then a certain limit h which is adjusted as a reconstruction parameter. Too 
large values of h lead to an unnecessary increase of CPU time since too many 
non- matching combinations are taken into account. On the contrary, too small 
values of h cause losses of correct combinations and consequently decrease the 
reconstruction efficiency. In our case, a reasonable value for h is 2 cm. 




Fig. 8. Crossing of planes representing two-dimensional solutions for three different 
orientations of wires. 

The determination of track parameters in three dimensions proceeds via min- 
imization of the function calculated on the basis of the differences between 
the measured distances of tracks from sense wires di and the distances of the 
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fitted tracks from the sense wires d: 



f. 



wliere tlie summation proceeds over all hits assigned to one track and 6di is 
an uncertainty of the measured distance di. In the case of the two-dimensional 
track reconstruction for wires oriented in one direction, the numerator in the 
above expression can be easily expressed in Cartesian coordinates using the 
formula: 



where a and b arc parameters in the straight line equation: x = az + b rep- 
resenting the track, x and z are coordinates perpendicular to the sense wires 
and {xi, Zi) is a point lying in the distance d^ from the sense wire in the lo- 
cation corresponding to the closest approach to the track (see Fig. 9). For 




Fig. 9. Definition of the point {xi, Zi) used in the two-dimensional track reconstruc- 
tion. 

the three-dimensional reconstruction, we use a coordinate system with the 
y-axis parallel to the vertical wires and the z-axis perpendicular to the detec- 
tion planes. The origin of the coordinate system is located in the first plane. 
Particle tracks are parametrized as straight lines in a conventional way: 



x'^a'z' + b', (3) 
y'^c'z' + d', (4) 

where a', b', c', d' are the searched track parameters. They can be linked to 
the parameters a, b corresponding to the two-dimensional solutions using the 
linear transformation: 
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a = a' COS ai + c' sin ai, (5) 
b — b' cos ai + d' sin ctj , (6) 

where ccj is the angle between the sense wires and the vertical direction. In- 
serting Eqs. (5) and (6) into (2) and the result into (1), the expression for 
takes the form: 

2 „ — (c;' COS r\ , + c'siiin,)-, — (//cosn, + c/' sin o;j)]^ 
^ " Y [1 + (a' cos ai + d sin Q;i)2] {bdif ' ^ ^ 



The derivative of the dimensionless term in the square bracket in the denom- 
inator is of the order of unity and is neghgible compared with corresponding 
derivative of the numerator divided by {ddi)"^ which is of the order of Zijbdi. 
Therefore, during the fitting procedure we regard the numerator just as a con- 
stant given by the initial values of the parameters and thus the minimization 
of is reduced to the hnear least squares problem. For the fitting we use 
the ,sudI/iHunction from "Numerical Recipes in C" [13]. The three-dimensional 
track fitting is performed for all combinations of two-dimensional track can- 
didates and the solution with the minimal value of is selected. The hits 
corresponding to this solution are then removed from the set of hits recorded 
in a given event and the reconstruction is repeated from the beginning until 
there are no further candidates of tracks. In this way multi-track events can 
be reconstructed. 



4 Calibration and reconstruction results 



The chamber is calibrated using the experimental data. In a first step an 
approximate drift time to drift distance relation is determined by integration 
of the drift time spectra as provided by the uniform irradiation method [12]. 
A typical spectrum obtained by summing up drift time spectra from all sense 
wires in one detection plane is shown in Fig. 10a. The width of the drift 
time distribution is about 300 ns which with the corresponding half of the 
cell width equal to 9 mm results in a mean drift velocity of 0.03 mm/ns. 
However, the space-time relationship obtained by integration of the drift time 
spectrum is not strictly linear (sec Fig. 10b); its slope increases slightly in the 
neighborhood of the sense wire. In the next step, corrections to this calibration 
are determined using an iterative procedure. For this, the average deviations 
between the measured and fitted distances of the tracks from the sense wires 
are calculated. These deviations are determined as a function of the drift time 
tj corresponding to the TDC channel j and the track entrance angle from 
the range (0°, 90°) divided into 9 intervals numerated by the index k: 

Ad{tj, 9k) =<d^{tj) > -d{tj, 9k), (8) 
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Fig. 10. Drift time spectrum for the first detection plane (a) and space-time re- 
lationship obtained by integration of the drift time spectrum in the time interval 
indicated by the arrows (b). The tail on the right hand side of the drift time spec- 
trum originates from electrons drifting from the space between the detection planes. 

where the average <> is taken over all hits which were recorded in the TDC 
channel j and correspond to the track entrance angle interval k. The space- 
time relation used in the first iteration is then corrected by the above deviation: 



After performing the reconstruction of tracks with the new space-time relation 
new corrections are calculated and so on. This procedure is repeated until 
the corrections become smaller than the position resolution of the chamber. 
This occurs in our case after 2-3 iterations. Fig. 11 shows differences Ad of 
the measured and fitted distances calculated as a function of the drift time 
for three subsequent iterations. The mean value of Ad deviates from zero 
only after the first iteration (upper panel in Fig. 11) and the corresponding 
correction to the space-time relation is of the order of 0.3 mm. For higher 
iterations the corrections are negligible. The standard deviation of Ad is about 
0.2 mm and is a measure of the single wire resolution. 

For testing the track reconstruction we used single tracks of protons scattered 
elastically on the COSY- 11 hydrogen target of the cluster-jet type at a beam 
momentum of 3.3 GeV/c. Events of the elastic proton-proton scattering were 
registered as coincidences between the scintillation hodoscope SI detecting 
the forward scattered protons and the monitor scintillator placed in the tar- 
get region on the opposite side of the beam, measuring the recoil protons. 
Additionally, we required firing of the seventh detection plane in order to as- 
sure that the proton tracks pass through the chamber. The detection efficiency 
per detection plane was estimated on the basis of multiplicities of coincidences 
between the detection planes. In 97% of cases all seven planes fire (see Fig. 
12) which means that the detection efficiency in a single detection plane is 



d'(tj,ek) = d{tj,ek) + Ad{tj,ek). 



(9) 
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Fig. 11. Differences of the measured and fitted distances calculated as a function 
of the drift time for the first detection plane and the angular bin 9 G (50° , 60° ) in 
three subsequent iterations of the calibration procedure. 

close to 100%. The multiplicity of hits shown in Fig. 13a is peaked around 
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Fig. 12. Multiplicity of coincidences between the detection planes for single tracks 
of elastically scattered protons. 

the value of 14 which is due to the fact that most of tracks are inclined with 
respect to the chamber and mostly two neighbouring cells from one detection 
plane fire for each track. For the reconstruction we used an "effective" position 
resolution of 220 fim for which the mean value of per degree of freedom is 
in the vicinity of 1. The reconstruction was treated as successful when the 
per degree of freedom was smaller than 5. This was the case for 95% of events 
triggered as elastic p — p scattering. For selected events in which all seven 
detection planes fired, the reconstruction was successful in 97%. In the recon- 
struction most of the registered hits were used for the fitting as one can see 
by comparison of Fig. 13a and Fig. 13b. The precision of the track parameters 
a, b, c, d was extracted from the diagonal elements of the corresponding covari- 
ance matrix. For typical track entrance angles of about 50° — 60° the resulting 
uncertainty of track position along the x- and ^/-axis is equal to 5h = 0.3 mm 
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Fig. 13. Multiplicity of hits for events of the elastic proton-proton scattering (a) 
and of hits preselected and used in successful reconstruction of single tracks (b) . 

and 5d — 3 mm, respectively. The uncertainty of the track inclination in the 
horizontal and vertical plane is about 1 mrad and 10 mrad, respectively. 



5 Conclusions 



A drift chamber with a c-shaped frame was constructed for the COSY-11 ex- 
perimental facility. The special shape of the frame allows for the detection 
of low momentum particles without disturbing the high momentum ejectiles 
which are registered in other detector components. This kind of chamber can 
be applied if no frame elements are allowed in the sensitive area of the detec- 
tion system. An example of a possible further application — besides the one 
discussed here — could be the detection of particles scattered at small angles 
with respect to the beam. In this case, two such chambers placed symmet- 
rically with respect to the beamline could be used for covering the forward 
angles. The chambers can be installed and removed without dismounting the 
beam-pipe contrary to an alternative solution of one chamber with a central 
hole for the beam-pipe. 

The track reconstruction algorithm which was developed for the chamber can 
also be used for other planar drift chambers containing cells with an electric 
field which has approximately cylindrical symmetry. It is required, however, 
that the chamber contains wires oriented at least in three different directions 
and for each direction there are at least two detection planes necessary. 

The chamber allows to determine the track position and inclination in the 
horizontal plane with an accuracy of about 0.3 mm and 1 mrad, respectively. 
In the vertical direction these accuracies are worse by about an order of mag- 
nitude, which is in accordance with the design values. For the tracking in the 
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vertical magnetic field of the COSY-11 dipole magnet no higher precision was 
envisaged. For other applications the precision can be improved by choosing 
a larger inclination of the wires. 
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